Gravitational Wave Astronomy Using Pulsars: 
Massive Black Hole Mergers & the Early Universe 

A White Paper for the Astronomy 8z Astrophysics Decadal Survey 



NANOGrav: 
The North American Nanohertz Observatory for 
Gravitational Waves 




Principal Authors: P. Demorest (NRAO, 434-244-6838, pdemores@nrao.edu); J. Lazio (NRL, 202-404- 
6329, Joseph.Lazio@nrl.navy.mil)] A. Lommen (Franklin & Marshall, 717-291-4136, andrea.lommen@fandm.edu) 
NANOGrav Members and Contributors: A. Archibald (McGill); Z. Arzoumanian (CRESST/USRA/NASA- 
GSFC); D. Backer (UC Berkeley); J. Cordes (Cornell); P. Demorest (NRAO); R. Ferdman (CNRS, France); 
P. Freire (NAIC); M. Gonzalez (UBC); R. Jenet (UTB/CGWA); V. Kaspi (McGill); V. Kondratiev (WVU); 
J. Lazio (NRL); A. Lommen (NANOGrav Chair, Franklin & Marshall); D. Lorimer (WVU); R. Lynch (Vir- 
ginia); M. McLaughhn (WVU); D. Nice (Bryn Mawr); S. Ransom (NRAO); R. Shannon (Cornell); X. Siemens 
(UW Milwaukee); I. Stairs (UBC); D. Stinebring (Oberlin) 

This white paper is endorsed by: ATA; LISA; NAIC; NRAO; SKA; US SKA; D. Reitze (LSC Spokesper- 
son, U Fl.); D. Shoemaker (LIGO Lab, MIT); S. Whitcomb (LIGO Lab, Caltech); R. Weiss (LIGO Lab, 
MIT) 



1 Science Opportunity: Exploring the Low-Frequency Gravita- 
tional Wave Spectrum 

Gravitational waves are fluctuations in the fabric of spacetime predicted by Einstein's theory 
of general relativity. Using a collection of millisecond pulsars as high-precision clocks, the 
nHz band of this radiation is likely to be detected within the next decade ( Jenet et al.||2005" ). 
The fundamental questions that will be addressed by these studies are: 

1. What is the nature of space and time? We suspect the local spacetime metric is 
perturbed by the cumulative effect of gravitational waves (GWs) emitted by numerous 
massive black hole (MBH) binaries. What is the energy density contained in this 
stochastic background of GWs? 

2. How did structure form in the Universe? Detection of GWs in the pulsar timing 
band will tell us whether MBHs formed through accretion and/or merger events. 

3. What is the structure of individual MBH binary systems? Recovering the 
gravitational waveform from individual systems will give us unprecedented insight. 

4. What contribution do cosmic strings make to the GW background (GWB)? 

The detection of cosmic strings would open a window into the early universe at a time 
inaccessible via the electromagnetic spectrum. 

5. What currently unknown sources of GW exist in the Universe? Every time 
a new piece of the electromagnetic spectrum has been opened up to observations (e.g. 
radio, X-rays, and 7-rays), new and entirely unexpected classes of objects have been 
discovered. 

The existence of GWs has already been inferred via the Nobel Prize- winning observations 



of the orbital decay of the PSR B1913+16 binary system (Hulse and Taylor 1975). While 
compelling and entirely consistent with general relativity, the behavior of this system offers 
only indirect evidence for GWs - the objective for the 21^* century is the direct detection 
and exploitation of GWs as a non-photonic probe of the Universe. 

Millisecond pulsars are old neutron stars that have been spun-up by mass accretion from 
a companion star to spin rates of hundreds of Hz. The rotational stability of these pulsars 
surpasses the majority of "normal" pulsars, and rivals that of atomic clocks. Pulsars emit 
a beam of radio waves that sweeps past the Earth once per rotation, appearing to us as 
a series of pulses. By precisely measuring the times of arrival of the the radio pulses on 
Earth, we can search for tiny perturbations due to GWs. This is distinct from detecting 
GWs emitted by the pulsars themselves. Rather, pulsar timing provides a means to detect 
any gravitational radiation crossing the Earth-pulsar line of sight, potentially from sources 
far outside our galaxy. 

All modern GW observatories operate on the principle that passing GWs cause tiny 
deviations in the distance between point masses. GWs are detected as changes in the light 
travel time between the points. In the case of the pulsar timing array (PTA) the path from 
each pulsar to Earth forms an arm of the GW detector. This detector is most sensitive to 
GWs with periods comparable to the total observation timespan, typically 1-10 years, which 
corresponds to nHz frequencies. The most desirable PTA involves millisecond pulsars evenly 
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distributed on the sky. A passing GW modifies the spacetime around the Earth in a manner 
that produces correlated shifts in the pulse times of arrival from the different pulsars. 

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav)[^ is 
an organization of astronomers, primarily from the U.S. and Canada, working to achieve 
GW detection using pulsar timing. Current projects include ongoing high-precision pulsar 
timing programs at Arecibo Observator}|^ and the Green Bank Telescope (GBT)[^ Here we 
describe the potential sources of low-frequency GWs and the current status and key 
advances needed for the detection and exploitation of GWs through pulsar timing (®. 



2 Science Context: Gravitational Wave Astrophysics 
2.1 Mergers of MBH Binaries 

It is now well established that mergers are an essential part of galaxy formation and evolution, 
and that massive black holes (MBHs, M > 10^ Mq) exist in the nuclei of most, if not all, 
large galaxies (see e.g. Ferrarese and Merritt 2000). Consequently, the product galaxy of 



many mergers will contain two MBHs. Due to dynamical friction, these two MBHs sink 
toward the center of the resulting galaxy's potential. 

As a MBH binary hardens, the strength of its GW emission increases. Once the system 
reaches a point where its semi-major axis is pc, GW emission becomes the dominant 
form of energy loss, and the two MBHs continue to spiral towards each other. Further, as 
the binary tightens, the frequency of the GW emission increases. For reference, an MBH 
binary with total mass M and semi-major axis a produces GWs with a frequency 
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Despite theoretical difficulties in fully understanding the production of these hard binary 
systems, observational results continue to provide evidence for their existence. One dramatic 
piece of evidence is the galaxy B0402-I-679, where high resolution radio imaging reveals two 



radio-loud nuclei separated by only 7 pc (Rodriguez et al. 2006). Follow-up observations of 



SDSS galaxies show strong evidence for MBH binaries with semi-major axes less than 1 kpc. 
A small number of these galaxies show [O ill] emission offset from the systemic redshift of the 
host galaxy. In two specific cases double-lined [O ill] profiles have been found, suggesting 
double AGN (Comerford et al. 2008). Most recently, a dual-broad-hne QSO system has 



been identified where the inferred MBH binary separation is only 0.1 pc (Boroson and Lauer 



2009). 



The general relativistic metric perturbation amplitude due to a GW is commonly pre- 
sented in terms of a dimensionless quantity called the characteristic strain he- The ensemble 
of MBH binaries is expected to produce a GWB whose amplitude spectrum has a power-law 
shape, hc{f) oc for a GW frequency /, where it is predicted that a = —2/3 (Phinney 
2001 Jaffe and Backer 2003). The strain amplitude is less certain, but expected to be in 
the range 10"^^ to 10"^^ at / = 1 yr'^ (I Jaffe and Backer! |2003l iSesana et al.||2008|). The 



^http : / / www . ricinograv . org 

^The Arecibo Observatory is a facility of the National Astronomy and Ionosphere Center, operated by 
Cornell University under a cooperative agreement with the National Science Foundation 

■^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc 
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Figure 1: Comparison of current and planned GW detectors, showing characteristic strain 
{he) sensitivity versus frequency along with expected source strengths. The Laser Interfer- 
ometer Gravitational Wave Observatory (LIGO), the Laser Interferometer Space Antenna 
(LISA) and PTAs occupy complementary parts of the GW spectrum. 



main contributors to this nHz GWB signal are MBH systems of mass M > 10^ Mq. Current 
pulsar timing experiments limit the GWB spectral amplitude to ^ 7 x 10~^^, depending 
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longer spans of data, PTA experiments become sensitive to lower GW frequencies, where the 
expected signal is stronger. The ongoing NANOGrav pulsar timing program will achieve a 
GW sensitivity well into the predicted he amplitude range in the next 3-5 years (Figures [l] 
andg. 

As these limits improve, and ultimately progress to a detection, they provide a new view 
of the history of MBH mergers throughout the Universe. Furthermore, measurements of 
the GWB spectral shape near 10 nHz could distinguish between various models of MBH 
binary formation (Sesana et al. 2008). LISA will be sensitive to the final MBH-MBH coa- 
lescence events for systems with M < 10^ M©. Pulsar timing arrays and LISA thus provide 
complementary views of these sources, encompassing the full range of MBH masses. 

2.2 Gravitational Wave Bursts and Individual Sources 

Pulsar timing has already been used to constrain GW emission from individual sources. For 
example, a proposed MBH binary within the radio galaxy 3C 66B (Sudou et al. 2003) was 
ruled out when its GW signature was not detected in existing millisecond pulsar data sets 
dJenet et al.]|2004D . 

GW bursts (events lasting less than a few years) are also potentially detectable. Sources 
might include highly eccentric systems near periapsis or the final inspiral of merging black 
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Figure 2: PTA sensitivity versus times for several scenarios: The current NANOGrav ob- 
serving program, and potential future PTAs of 3 and 10 times better GW sensitivity. The 



shaded area shows the expected amplitude range of the MBH-MBH GW signal (see ^2.1). 



holes (Lommen et al. 2009 ). GW waveforms of such events would encode detailed information 



about the burst source, such as the masses and spins of the inspiraling black holes. 
2.3 Cosmic Strings and Exotica 

Pulsar timing experiments may provide a unique window into particle physics at the high- 
est energy scales. Cosmic strings, theorized line-like topological defects, may form during 
phase transitions in the early Universe, due to the rapid cooling that took place after the 



Big Bang (Kibble 1976). Recently it was shown that cosmic string production is generic in 



supersymmetric grand unified theories (Jeannerot et al. 2003). Furthermore, string theoret 



ical cosmology predicts the formation of so-called cosmic superstrings, different from regular 



field theoretical cosmic strings (Polchinski 2005). 



Cosmic strings and superstrings are expected to produce a stochastic GWB analogous 



to the cosmic microwave background, as well as bursts of GWs (Damour and Vilenkin 2001 



2005 Siemens et al.||2006 ). Because of their sensitivity at very low frequencies, pulsar timing 
arrays place the best constraints on viable cosmic string models (Siemens et al. 2007). With 
increased sensitivity, PTAs could detect cosmic (super) strings. 

Pulsar timing experiments could result in the detection of other exotica. Gravitational 
waves are a means for probing the fundamental structure of the space-time of the Universe. 
With less than 20% of the matter in the Universe emitting electromagnetic radiation, we are 
likely to be surprised by what we "see" in GWs, the generation of which is caused directly 
by the movement of mass, not the coupling to the electromagnetic force. 

3 Key Advances For A Pulsar Gravitational Wave Observatory 

The sensitivity of a PTA is determined by the number and distribution of the pulsars under 
observation, the cadence with which they are observed, and the precision with which the 
pulse times of arrival are measured. 

Pulsar timing precision is quantified by the root-mean-square (RMS) residual pulse arrival 
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Figure 3: Published RMS pulsar timing residuals versus time, showing exponential improve- 
ment (Demorest & Jenet, 2009) which positions us to detect GWs within the next decade. 

time after a fit to a standard model of pulsar rotation, binary motion. Earth motion, and 
interstellar propagation effects. The RMS residual for a given source is determined by its 
flux density, characteristic pulse shape, emission stability, rotation stability, the scintillation 
and scattering of its signal as it traverses the interstellar medium, and the radio telescope 
equipment used to observe it (telescope area, system temperature Tgys, bandwidth, data 
acquisition instrumentation, and detection algorithms). Currently there are several pulsars 
with RMS residuals approaching 100 ns, and roughly 20 more with residuals less than 1 /is. 



Jenet et al. (2005) showed that with 100-ns level timing on 20 pulsars, the stochastic GWB 
is detectable in ~5 years. Realizing this goal will come from a two-pronged approach: We 
must flnd additional pulsars suitable for high-precision timing and also improve the timing 
precision of known sources. 



Demorest and Jenet (2009) have recently suggested that published RMS pulsar timing 



residuals over the past two decades show an exponential improvement with time analogous 
to Moore's Law for computer processors (Figure [s]), improving by a factor of 2 every ~3 
years. Several observational advances are required if we are to continue this trend: 
A. Pulsar Surveys The sensitivity of a PTA scales directly with the number of pulsars 



in the array (Jenet et al. 2005). Given sufficient observing time, we expect to be make a 



detection of GWs with currently known pulsars within 5 years. However, to fully characterize 
the gravitational waveforms and to maximize the scientific return, more high-timing precision 
millisecond pulsars are needed. It is especially important to include more pulsars in directions 
widely separated from the current set of objects. Three current searches are now turning 
up such objects: the PALFA L-band multibeam survey at Arecibo, the GET 350 MHz drift 
scan survey, and the new Parkes L-band Digital Survey. In the next decade, a new GET low- 
frequency pulsar survey would be particularly advantageous because it would increase the 
number of Northern hemisphere pulsars. This area of the sky is currently under- represented 
in timing arrays. As pulsar flux increases at low radio frequencies (from ~1 mJy at lA GHz 
to ~10 mJy at 400 MHz), this survey will identify many new pulsars. 
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Table 1: International PTA telescope time in terms of a 100-m dish with T, 
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^ Includes the effects of reflector efRciency and partial illumination. 

^ This represents the combined observing time of four European 100-m class dishes. 

Equivalent single-dish diameter. 

20 pulsars with <100 ns RMS timing. 
" >40 pulsars with ;^100 ns RMS timing. 

B. Time on High Sensitivity Facilities To achieve the necessary precision on these 
faint objects we require multifrequency observations on 100-m class or larger radio tele- 
scopes. Four-frequency observations are necessary in order to precisely fit for the frequency- 
dependent dispersion caused by free electrons in the ISM, which typically produces per- 
turbations of ~/is. Table [T] illustrates the total worldwide telescope time being devoted to 
pulsar timing with the goal of GW detection, weighted by telescope sensitivity. For refer- 
ence, several future facilities are listed: By 2011 five European telescopes will be combined 
into a phased array, the Large European Array for Pulsars (LEAP), with an effective area 
equivalent to a 200-m single dish. Other future facilities are the full 350-dish Allen Telescope 
Array (ATA-350), the Expanded Very Large Array (EVLA) and the Square Kilometer Array 
(SKA). 

The table shows that world resources currently provide about 300 100-m hours per month. 
We estimate that at minimum, GW detection requires ~500 100-m hours per month, based 
on observing 20 pulsars every 2 weeks, for 3 hours at each of 4 radio frequencies in order 
to obtain 100-ns or better timing precision. To fully characterize the GW sources requires 
at least twice as many pulsars (Lee et al. 2008), and an effort to upgrade receivers and 



backend instrumentation to handle ~1 GHz total bandwidth. This process is underway at 
many existing radio telescopes. Clearly, to fulfill the goal of detecting and characterizing the 
stochastic GWB, as well as continuous and burst sources, we require more resources than 
are currently available. The table suggests that the ATA and EVLA could possibly provide 
the additional time needed in the near future, and the SKA farther in the future. 

NANOGrav, the Parkes Pulsar Timing Array (PPTA) and the EPTA are in the process 
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of organizing themselves as the International Pulsar Timing Array (IPTA) for the purpose of 
optimizing these international resources. In addition to sensitivity, a full optimiziation must 
consider such additional factors as telescope sky coverage, frequency agility, and backend 
instrumentation. 

C. Algorithm Development Algorithms for both radio pulse detection and GW detection 
must be improved. With increased telescope sensitivity, pulse arrival times are increasingly 
susceptible to systematic effects. Among the algorithms to be developed are methods for 



characterizing and compensating for the effects of the interstellar medium (Foster and Backer 



1990 



et al. 



You et aL]|2007); methods for effectively mitigating radio frequency interference (Stairs 



2000); and methods for fully utiziling the available polarization information (van 



Straten 2006). Algorithms for optimal extraction of the GW signal will be built upon the 



recent advances of LIGO and LISA data analysis; some of this work has already begun ( Jenet 



et al. 2005; van Haasteren et al. 2008; Anholm et al. 2008 Lommen et al. 2009). 



4 Summary 

Given sufficient resources, we expect to detect GWs through the IPTA within the next five 
years. We also expect to gain new astrophysical insight on the detected sources and, for the 
first time, characterize the universe in this completely new regime. The international effort 
is well on its way to achieving its goals. With sustained effort, and sufficient resources, this 
work is poised to offer a new window into the Universe by 2020. 
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